Charge pairing and superconductivity in high-Tc cuprate superconductors 



Eduardo C. Marino and Marcello B. Silva Neto 
Instituto de Fisica, Universidade Federal do Rio de Janeiro, Caixa Postal 68528, Rio de Janeiro 

(February 1, 2008) 



RJ, 21945-970, Brazil 



O 
O 

(N 

> 
O 

o 



I 

o 
o 



> 



o 



o 
o 



X 



We propose a model for high-Tc superconductors that includes both the spin fluctuations of the 
Cu"'"'' magnetic ions and of the spins of O doped holes (spinous). The charge of the dopants 
(holons) is associated to quantum skyrmion excitations of the Cu++ spin background. The quantum 
skyrmion effective interaction potential is evaluated as a function of doping and temperature, indi- 
cating that Cooper pair formation is determined by the competition between these two types of spin 
fluctuations. The superconducting transition occurs when the effective potential allows for skyrmion 
bound states. Our theoretical predictions for the superconducting phase diagram of La2-a;Srj;Cu04 
and YBa2Cu306+i: are in good agreement with experiment. 



PACS number(s): 74.72.Bk, 74.25. Ha 

Introduction. High-temperature underdoped cuprates 
exhibit a wide variety of interesting physical phenomena, 
like Neel and metal-insulator transitions, non- Fermi liq- 
uid behavior, pseudogap, etc., and have inspired a large 
amount of theoretical and experimental work for about 
fifteen years. In spite of that, even the nature of the 
ground state and of its elementary excitations, have not 
yet been fully determined and many different pictures are 
availabe 

Another fundamental point yet to be understood is 
the mechanism of charge pairing. It is by now well es- 
tablished that antiferromagnetic spin correlations play an 
important role in the dynamics of the system, even after 
the destruction of the Neel state. Indeed, different spin- 
fluctuation models have been successfully used to explain 
the observed spectral weight in ARPES data of high-Tc 
materials [Q, as well as other anomalies Moreover, 
the idea of spin-fluctuation induced charge pairing and 
superconductivity has been used recurrently [Q. 

In this work we propose a theory for high-Tc cuprates 
that takes into account the spin fluctuations of the Cu++ 
magnetic ions and of the O doped holes as independent 
degrees of freedom. The charge of the dopants (holons) 
is associated to skyrmion quantum spin excitations of 
the Cu"^"*" background, which in the Neel phase are fi- 
nite energy defects closely related to their classic coun- 
terparts whereas in the quantum disordered phase are 
nontrivial zero energy purely quantum mechanical exci- 
tations. The spin of the doped holes (spinous), on the 
other hand, is represented by chargeless, massless Dirac 
fermion fields We calculate the effective interaction 
potential between these quantum skyrmion topological 
excitations in order to study charge pairing. It becomes 
clear that Cooper pairing is controlled by the competi- 
tion between the spin fluctuations of Cu^"*" magnetic ions 
and those of the doped holes. Our predictions for 

the Tsc line are in good agreement with experiment for 
both La2-a;Sra:Cu04 and YBa2Cu306-i-a;. 



The model. Our starting point will be the generalized 
spin-fermion model described by the Hamiltonian 

{i,j),a i,a 

+ JJ2S.-S,, (1) 

which arises from the strong coupling limit of the three 
band Hubbard Model (3BHM) §]. In the above expres- 
sion. Si represent the localized spins of Copper ions, 
which interact through the superexchange J, ^, a = 
l..A^ = 2, is the hole creation operator, tp is the hopping 
term for holes, Jk is a Kondo like coupling between the 
spins of Cu"'""'" ions and the spins of O holes, and we 
have retained the usually ignored onsite Coulomb repul- 
sion between O holes. Up with ni.Q, = cl^Ci,a- 
The reason is that realistic estimates from the 3BHM 
suggest that Up/tp ~ 10 0, being rather large, and thus 
we can perform a tp/Up expansion. Second order pertur- 
bation theory in tp/Up will give rise to a superexchange 
Jp = 2tp/Up between oxygen spins and we end up with 
a. t — J model for the holes. 

The mean field (large N) solutions of the < — J model 
are well known and it has been established that a tt- 
flux phase has minimum energy, at least at the saddle- 
point level {N oo) We can write the electron 
in terms of a charged spinless boson fii (holon) and a 
chargeless spin-1/2 fermion fi^a (spinon), and, as usual, 
we decouple the four particle interactions by introduc- 
ing the d-wave auxiliary fields Xij = (^fiafj,a^ ^^'^ 

^ij = {fi,]fj,l - fidfj,])^ which are nonzero for T < T*, 
where T* is the pseudogap temperature. If we then 

neglect charge fluctuations, ^^^^j^ ~ = const., 

we find that the lowest lying excitations of the 7r-flux 
phase are massless, chargeless, spin carrying Dirac Fermi 
fields 1^ whose dynamics is described by the Lagrangian 
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vpj ■ V] 4'a,\, where A = 1,2 



label the two Fermi points at (7r/2, ±7r/2), ~ {dr, V), 
Ifi = (70,7) = (io-^,o-a;,cry), w_F = 2ax is the dopant 
Fermi velocity (with a being the lattice spacing and x 

Ja,\. 

for (o)dd and (e)ven lattice sites. The long wave- 
length fluctuations of the localized Cu"'"'" spins, on the 
other hand, are described by the CP^~^ Lagrangian 



the constant amplitude of \xij\) ipa.x 
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CP'' 



with z. 



(l/2go) |(i9m - ^■^t^)zi\ , where S = z, 
, i = 1..N 



fields such that zjzi 



— 2, being Schwinger boson 
1, Af_, = -izid^Zi, and go is 
a bare coupling constant. It is now convenient to per- 
form the local canonical transformation f/' Utp, where 

Zl -Z2 



U = exp 9 (^"^ ^ '^^(^)' ^'^d q is arbitrary. 

Now the Kondo coupling term in (|l|) reduces to a chem- 
ical potential term, since WS ■ all = Uz. Also, since 
U^dfJJ — iqazA^+ negligible nondiagonal terms, we end 
up with the effective theory 



Z = / VzVzViiVipVAf, 5[zz - l]e 



-s 



(2) 



where 

5* = 



1(9^ - iAf,)zi\ 



E 

Q = l..Af,A=l,2 



(3) 



Holons as quantum skyrmions. In previous works 
we have proposed a model for doping quantum Heisen- 
berg antiferromagnets, that successfully described the 
magnetization curves and the AF part of the phase dia- 
grams of both LSCO and YBCO. One of the important 
consequences of that model was the observation that each 
hole added to the Cu02 planes creates a skyrmion topo- 
logical defect on the Cu++ spin background, in agree- 
ment with earlier proposals The dopant charge, in 
particular, was found to be attached to the skyrmion 
charge and consequently its dynamics becomes totally 
determined by the quantum skyrmion correlation func- 
tions. Despite the fact that the model proposed in Q is 
restricted to the antiferromagnetic part of the phase di- 
agram, we shall nevertheless pursue the picture in which 
skyrmions are in general the charge carriers of the doped 
holes. This will allow us to treat the bosonic variable 
introduced above as a quantum skyrmion operator. In 
particular, we shall exploit this idea in the quantum dis- 
ordered phase, S > Saf, where the skyrmions are purely 
quantum mechanical and have zero energy. 

The full treatment of the quantum skyrmions of the 
theory described by (0) has been carried out in In 



the renormalized classical regime, go < gc {gc = Stt/A), 
we have 



g-27rpja;-i/| 

|x-y|9V2 ' 



(4) 



where ps — 1/ go — 1/gc- Conversely, for the theory stud- 
ied in j|] the corresponding correlator was found to be 



(5) 



where the expressions for ps{6) and a(6) have been 
carefully determined in [p[. In particular, a{S) — 
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{n6)^ with n = 1 for YBCO and n = 4 
for LSCO, the factor of four being a consequence of the 
existence of four branches in the Fermi surface for this 
compound, as discussed in a^M is the electromag- 
netic fine structure constant and we see that the contri- 
bution from the electromagnetic coupling is negligible. 
The Ps{S) function is given by Ps{S) = ps(0)[l — AS"^], for 
YBCO and psiS) = Ps(0)[l - B6 - CS"^]^/"^, for LSCO, 
and again the different behavior being ascribed to the 
form of the Fermi surface in each case The con- 
stants A^B and C have been evaluated from first prin- 
ciples in 1^. In order to obtain the ^-dependence of 
the spin stiffness ps and of the spinon coupling q in our 
model (|), we now match the two correlation functions 
in (li) and (ordered phase), obtaining ps = Ps{S) and 



q = 



r 128 - 



1/2 



(nS). The sublattice magnetization in the 

ordered phase is given by M{5) — ^ p{S), and conse- 
quently Saf can be obtained from p{Saf) — 0, both in 
good agreement with experiment, see For S > Saf, 
on the other hand, where ps — 0, we shall assume that 
the expression for q{6) still holds. This is quite reason- 
able since q was introduced by a local canonical trans- 
formation, and at least locally there is short range AF 
order. 

Cooper pair formation. Let us now investigate the 
conditions for Cooper pairing. We shall first intro- 
duce in the partition function (p|) the skyrmion current, 



Jf" = ^e^"^aa^;3, through the identity 



(6) 



Integrating over zl,Zi and ?/'a:'0a, we obtain, at leading 
order, the effective Lagrangian 

N 

^e//[-4^] = Y^^(x,T)n^"^(x-y,T-r')A(y,r'), (7) 

where 11'^'^ (x — y, r — r') has Fourier transform given 
by n'"''(p,ie™) = n^''(p,ie„0 -f n^''(p, ie„j), which are 
respectively the contributions to the finite temperature 
vacuum polarization coming from the complex scalar 
fields Zi (Schwinger bosons) and fermions V'a,A (spinous). 
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In order to obtain the effective current-current interac- 
tion between skyrmions, we use an exponential represen- 
tation for the i5-function in and integrate over Afj, and 
the corresponding Lagrange multipUer field. The result 
is 



Z 



where = H^^" (p) / p'^ , x — (r, x) and p — (ie„j,p). 

The real time effective interaction energy between static 
skyrmions (e™ — 0) is then 

Hi = 27r2 j d^x j A^y p(x) E""(x - y; 0) p(y), (9) 

where p(x) = j7o(x) is the dopant charge density and 
I]°'^(x — y;0) has Fourier transform given by I]''°(p) — 
nB(p) + nF(p) with 



A 



1 



27r 27r 



da; v^|ppa;(l — x) + rri^ 



X coth 



■y/|ppa;(l — x) + m? 
2fc^T 



and 



I{f{-p) = ^ dx ^/\p\^x{l-x) tanh ^ 



(10) 



Vlpp^(i-^) 



2fcsT 



(11) 



In the above expressions, m is the inverse correlation 
length of the quantum disordered phase of the CP^^^ 
model, where A = Stt (1/gc — I/50) and ps = 0. At order 
N, it is given exactly by m{T) A-K2fcsTe-^/'=«^ [0. 

For two charges at positions xi and X2, we have 
p(x) = (5(^'(x — xi) -I- (5^^^ (x — X2). After discarding self- 
interactions, we obtain (r = Xi — X2) 



V{r) = / d2pS""(p) e'P-'' + Vi{r) 



(12) 



where we have also introduced the centrifugal barrier po- 
tential between the two charges that form the Cooper 
pair, Vi{r) = l{l + l)n'^ /2M*r'^ , with I specifying the rel- 
ative orbital angular momentum of the pair and M* the 
effective mass of the charges. 

Zero temperature limit. It is well known that in high- 
Tc cuprates, Cooper pairs form at relatively short dis- 
tances. In this limit, large |p|, we have, at T = 



V{r) 





■ 1 2q^' 







+ Vi{r). (13) 



The above expression clearly shows a competition be- 
tween the spin fluctuations of the Cu^"*" spins (first term) 
and of the doped spins (second term). For small 



enough doping, q"^ < 1/2, the potential is always repul- 
sive and there is no charge pairing. For > 1/2, on 
the other hand, the potential has a minimum and charge 
(skyrmion) pairing occurs. We conclude that the critical 
doping for the onset of superconductivity is determined 
by the condition q^(^sc) = 1/2. We observe that with- 
out the Cu^^ background, the interaction potential ( p^ ) 
would always have bound states for any g ^ 0, at zero 
temperature, and Ssc — 0- This is what happens in the 
mean field phase diagram of Kothar and Liu (l|]. We 
see that the effect of the Cu"*""*" background is to shift the 
value of Ssc to its correct position in the phase diagram. 

Determination of Ssc- From the expression of q in 
terms of 5 {q = ^/2a{6)), we see that 6sc is only deter- 
mined by the shape of the Fermi surface of the compound. 
In particular, we see that S^c'^'^ = 4(55^'^*^, a result that 
is verified by experiments, if we take in account the rela- 
tion between S and the stoichiometric doping parameter 
X, namely S = x for LSCO and S = x - 0.20 for YBCO. 
Another prediction of our model is that compounds with 
similar Fermi surfaces should have the same supercon- 
ducting critical doping Ssc- We get S^^^'^ = 0.318 and 
SscF'^ = 0.079, which have a fairly good agreement with 
experiment. 

Disorder. Disorder may be modelled in the ordered 
Neel phase of a doped antiferromagnet by considering a 
continuous random distribution of spin stiffnesses [ p^ . 
The effect of introducing a Gaussian x|p — Psl""^ distri- 
bution, with exponentially suppressed magnetic dilution, 
in the original model is a correction for a{5), namely 
a{5) a' {5) = a[5) + v Choosing u = \ for 

both compounds, we get Ssc = ^ \J ^^^2^ 1 or equivalently 

x^^'^o = 0.425 and x^^P^ = 0.056, in good agreement 
with experiment. 

Finite temperature. For lis we shall expand in 
kBT/m, since it is clear that m{T) > kBT,^T. Hb will 
then be simply given by its zero temperature limit, where 
TO = A. For lip, on the other hand, such a low T ex- 
pansion is not necessarily valid even for |p| ^ UbT. We 
will then have to split the integral over the Feynman pa- 
rameter x in ( pT| ) into three parts. For Q < x < Xc and 
1 — Xc < X < 1, Xc — (fcsT/lpl)^, we will use a high T 
expansion, while for < a; < 1 — we use the low T 
expression. We obtain 



EOO(p)=(l-^'^') 



IQq^T'^ - 4to3 



8|P| 



7r|p|- 



2|P|^ 



3^|p|' 



(14) 



Inserting this in (p^), we get ^(r), and from the tresh- 
old conditions for the formation of bound states, namely 
V'{yq) = and V"{yq) = 0, we obtain the relation 



(kBTsc) 



3 7r(l — 2q^)a^ ma^ 3Trm^a 



512^2 



32g2 12852 4g2 
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where a — Hvp/tq, m — A, and Tq is the mmimum of the 
potential (it also measures the size of the Cooper pair). 
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FIG. 1. Plot of curve ( jl^ ) for La2-a;Sri,Cu04. Experimen- 
tal data from l|l3l . 



new saddle-point. In the second plot (YBCO) we have 
shifted the curve (dashed part) to the right in the regions 
S = [0.52, 0.7] and 5 = [0.8, 1] in order to comply with the 
effects of the out-of-plane 0-Cu-O chains, which produce 
the observed 60 K and 90 K plateaus, where the extra 
holes do not enter in the Cu02 planes. Furthermore, for 
YBCO, T* is higher than T^ax (S = 1) and therefore 
imposes no restrictions to our results. 

Final remarks. We would like to remark that our the- 
ory (|^) also gives a simple interpretation for the pseu- 
dogap phenomena. Indeed, for Tsc < T < T* spinous 
are paired into d-wave singlets but holons (skyrmions) 
repel each other and there is no superconducting state. 
Only for T < Tsc we do have Cooper pair formation and 
superconductivity. 
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FIG. 2. Plot of curve ([ij) for YBa2Cu306+^. Experimen- 
tal data from 



Comparison with experiment. In order to make contact 
with experimental data we need the doping dependence 
of A. For YBCO, we use A{d) = AqUS/Saf)^ - 1], in 
agreement with the results of [^, with Aq = 8.8 meV. 
For LSCO, we shall use an expression that fits the exper- 
imental data of [0, namely A((5) = Aa[{S/SAF)^ -l]^^^ , 
with Ao ~ 1.0 meV. For the T — AF quantum critical 
point 6af, we know from experiments that Saf = 0.22 
for YBCO and 6af = 0.02 for LSCO. Inserting in (|l|) 
the values of 6sc a.t T — 0, obtained previously, we get 
a relation that fixes M*Vp with respect to tq. In figs. 
|l| and I we plot the curve (jTs]) for LSCO and YBCO, 
respectively, with tq = 38 A, Kvf = 0.18 eV A for LSCO 
and hvF = 1.08 eV A for YBCO, and Z = 2 (d-vave 
pairing). In the first plot (LSCO), the dashed part is 
in the region where T > T* and we should move to a 
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